Evidence of Urban-Induced Precipitation Variability in 
Arid Climate Regimes 

J. Marshall Shepherd 
NASA/GSFC 


Submitted to Journal of Arid Environments 
May 2005 


Corresponding Author: 

Dr. J. Marshall Shepherd 

Official and Permanent Address: 

NASA Goddard Space Flight Center 

Code 613.3 

Greenbelt, MD 20771 

Email: marshal 1 . shepherd@nasa. gov 

V: (301) 614-6327 

F: (301)614-5492 


Abstract 

The study employs a 108-year precipitation historical data record, global climate observing network 
observations and satellite data to identify possible anomalies in rainfall in and around two major arid 
urban areas. Phoenix, Arizona and Riyadh, Saudi Arabia. The analysis reveals that during the monsoon 
season, locations in northeastern suburbs and exurbs of the Phoenix metropolitan area have experienced 
statistically significant increases in mean precipitation of 12 to 14 percent from a pre-urban (1895-1949) 
to post-urban (1950-2003) period. Further analysis of satellite-based rainfall rates suggests the existence 
of the anomaly region over a seven-year period. The anomaly cannot simply be attributed to maximum 
topographic relief and is hypothesized to be related to urban-topographic interactions and possibly 
irrigation moisture. Temperature records suggest that Riyadh has experienced an adjustment in mean 
temperature in response to the growth of urban surfaces (e.g., the so-called urban heat island effect). 
While ground-based precipitation records also indicate an upward trend in mean and total precipitation in 
and around Riyadh in the last 10-15 years, it is difficult to attribute the increase to urbanization because 
other less urbanized stations in Saudi Arabia also show a similar increase. Recent satellite-based 
precipitation estimates indicate an anomaly region 50-100 km north of Riyadh, but this study is not robust 
enough to conclusively link it to urbanization although certain climate-regime attributes suggests that it 
might be. 
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I. Introduction 


Water is essential to life in the Earth system. The water cycle components that sustain life are becoming more 
scarce and polluted. The most recent (1999-2004) drought experienced in the southwestern United States is the 
seventh worst in the approximately 500-year proxy tree-ring record (Piechota et al. 2004). As a result, many 
regions contemplated “drought emergencies” in which severe water restrictions are implemented. Though large- 
scale forcing likely controls drought processes (Hidalgo 2004), there is increasing evidence that anthropogenic or 
“human-related” activities can significantly alter precipitation processes. Urbanization is an example of 
anthropogenic forcing. Recent studies (Inoue and Kimura 2004; Burian and Shepherd 2004; Molders and Olson 
2004; Dixon and Mote 2003; Shepherd and Burian; Bomstein and Lin 2000) continue to provide evidence that urban 
environments can modify or induce precipitation under a specific set of conditions. Results of these recent studies 
are consistent with previous historical work in this research area (Changnon et al. 1981; Huff 1986). 

Arid and semi-arid regions of the southwestern United States and other parts of the world are rapidly 
developing and placing greater demands on the environmental system. In the past fifty years, Phoenix has expanded 
from a predominantly agricultural center to an urbanized region with extent 700 percent larger than its size in the 
middle of the twentieth century (Rex 2000). Riyadh’s population grew from about a half million people in 1972 to 
almost two million by 2000. Saudi Arabia experienced urbanization later than many other countries; in the early 
1970s its urban-rural ratio was still about 1:3. By 1990 the ratio had reversed to about 3:1. In the mid-1970s 
Riyadh’s population was increasing by about 10 percent a year. Irrigation also significantly increased between 1972 
and 1990 (Metz 1993) southeast of Riyadh. 

Figure 1 is a LANDSAT image depicting the rapid growth of Phoenix over the past few decades. Figure 2 is a 
combined LANDSAT and ASTER (Advanced Spaceborne Thennal Emission and Reflection Radiometer) image 
depicting the rapid urban growth of Riyadh, Saudi Arabia. ASTER is an imaging instrument flying on Terra, a 
satellite launched in December 1999 as part of NASA’s Earth Observing System (EOS). There are a sampling of 
studies that discuss the anthropogenic impacts on the environment in arid regimes such as Kuwait City, Kuwait 
(Akber et al. 2001), Riyadh, Saudi Arabia (e.g., El-Sharif 1985; Modaihsh 1997), and Tehran, Iran (Karamouz et al. 
2001). Additional localized impacts resulting from urban and agricultural development have been documented, 
including the modification of local climates (e.g., Qurcshi and Khan 1994). 


However, there have been relatively few studies that focus on arid urban areas and precipitation enhancement. 
Balling and Brazel (1987) observed more frequent late afternoon storms in Phoenix during recent years of explosive 
population growth. Diem and Brown (2003) (hereafter DB03) suggested that anthropogenic activities in the Phoenix 
area affect summer precipitation totals, particularly the Lower Verde basin zone. Because arid regions lack 
substantial rainfall, investigators have not focused on these regimes for study of precipitation variability linked to 
human activity. Yet, they offer very unique attributes that are ideal for studying the problem. First, many areas in 
arid regimes have undergone rapid growth in the last 20-30 years which enables assessment of pre- and post- urban 
growth and precipitation data using a wealth of robust ground and satellite data sources. Second, the aridity of 
Phoenix and Riyadh minimizes the risk of precipitation contaminants from large-scale weather forcing such as fronts 
or low pressure systems thereby enabling a focus on convective-regional scale process linked to urbanization. Third, 
arid cities rely heavily on irrigation to support human activity. Increased irrigation could lead to moister 
environments to support or enhance precipitation processes. 

2. Science Objectives, Materials, and Methodology 

a. Phoenix Study 

We chose the cities of Phoenix, Arizona and Riyadh, Saudia Arabia for the study for the aforementioned 
reasons. The methodology for each city differs because of data availability. DB03 offered very compelling findings 
but their study only investigated the relative increases in precipitation totals for a sample of stations in the basins 
east of Phoenix and only utilized precipitation data from 1950 to 2000, which is essentially the post-urban period. 
Herein, urban processes are defined as including urban land use, urban aerosols, and human irrigation procedures. 
Our work extends DB03 in several important ways. 

First, we utilize a high-quality, topographically-sensitive, 108-year precipitation data record provided by the 
Spatial Climate Analysis Service (SCAS) at Oregon State University. This 4-km resolution dataset resolves the pre- 
urban (1895-1949) and post-urban precipitation climatology for the southwestern U.S. This capability is necessary 
for detecting pre-urban and post-urban effects since the city cannot be physical removed for experimentation. 
PRISM (Parameter-Regression Independent Slopes Model) (Daly et al. 2002) is used to spatially distribute ground- 
based observations. The primary sources of observations currently used by SCAS for mapping precipitation are: 1) 
COOP (NWS COOPerative observations), 2) SNOTEL (USDA NRCS SNOw TELemetry), and 3) ASOS 
(NWS/FAA Automated Surface Observing System). The total number of precipitation stations was approximately 


12600 (figure 2). A particularly attractive feature of PRISM for studying mountainous regions is its ability to 
calculate linear climate-elevation relationships. 

Second, we investigate the relative magnitude of mean precipitation totals in Phoenix and immediate 
surrounding regions. The goal is to verify results presented by DB03 that indicated a tendency for elevated 
precipitation totals in the Lower Verde basin region. We extend DB03 by integrating analysis of stations in the 
western and northern surrounding area as well as the Phoenix urban area. Statistical significance tests were 
conducted for the pre- and post-urban analyses and the relative magnitude analysis. 

The experimental Tropical Rainfall Measuring Mission (TRMM) 3B42 real-time rainfall product (TRP, version 
5) is used to examine mean rainfall rates of the central Arizona region during the extreme drought conditions of 
2003. The goal is to detect the precipitation anomaly under relatively homogeneous drought conditions. TRP 
provides rainfall rate and accumulation estimates at 0.25°-resolution every three hours over the latitudinal range of 
60°N - 60°S. TRP is a merger of all available SSM/I and TRMM microwave imager precipitation retrievals into a 
"high-quality" (HQ) precipitation estimate. The product also included precipitation estimates from geostationary 
infrared observations that are calibrated by the HP precipitation estimate (Huffman et al. 2003). Satellite data is 
useful because quantitative analysis of WSR-88D radar data is limited in terrain-intensive regions around Phoenix. 

b. Riyadh Study 

For this component of the study, we use data from the Global Climate Observing System (GCOS) Surface 
Network (GSN) and previously published data from Alkolibi (2002). The Global Climate Observing System 
(GCOS) was established in 1992 to ensure that the observations and information needed to address climate-related 
issues are obtained and made available to all potential users. It is co-sponsored by the World Meteorological 
Organization (WMO), the Intergovernmental Oceanographic Commission ( IOC ) of UNESCO, the United Nations 
Environment Programme ( UNEP ) and the International Council for Science ( ICSU ). GCOS is intended to be a long- 
term, user-driven operational system capable of providing the comprehensive observations required for monitoring 
the climate system, for detecting and attributing climate change, for assessing the impacts of climate variability and 
change, and for supporting research toward improved understanding, modeling and prediction of the climate system. 
It addresses the total climate system including physical, chemical and biological properties, and atmospheric, 
oceanic, hydrologic, cryospheric and terrestrial processes. The National Oceanic and Atmospheric Administration's 
(NOAA) National Climatic Data Center (NCDC) is responsible for building a permanent data base of historical 
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GSN daily and monthly data submissions, along with the appropriate station metadata history, and for providing free 
and open user access to this information via the World Wide Web. Our access point for the data was 
http ://ww w . ncdc . noaa. gov/oa/c 1 imate/gsn/ gsn i n fo.htm 1 . The study also employs data from version 6 of the TRMM 
3B43 algorithm. The purpose of Algorithm 3B-43 is to produce the "Tropical Rainfall Measuring Mission (TRMM) 
and Other Data" best-estimate precipitation rate and root-mean-square (RMS) precipitation-error estimates. These 
gridded estimates are on a calendar month temporal resolution and a 0.25-degree by 0.25-degree spatial resolution 
global band extending from 50 degrees South to 50 degrees North latitude. The global rainfall algorithm (3B-43) 
combines the estimates generated by combined instrument rain calibration (3B-42) and global gridded rain gauge 
data from NOAA’s Climate Prediction Center and/or global rain gauge product, produced by the Global 
Precipitation Climatology Center (GPCC). For more details of the algorithm, go to 
http trm m . gsfc . n asa , gov/3 h4 3.html. 

The approach is to identify, using historical gauge data, any trends in temperature and precipitation in and 
around Riyadh and other Saudi Arabian stations in the GSN. Then, using the TRMM-based product, identify 
possibly precipitation anomalies that may be linked to the urbanized Riyadh environment. 

c. Climatology and Topography 
i. Phoenix 

The southwest monsoon influence is greatest over portions of Mexico and extends into southeastern Arizona 
where average precipitation during July, August, and September exceeds 50% of the mean annual total. The Phoenix 
metropolitan area receives approximately 40% of its mean annual rainfall during the monsoon period. During the 
monsoon period, near-surface moisture availability in the Phoenix region peaks. The mean Tucson sounding 
(Wallace et al. 1999) in figure 3 indicates that on active monsoon days, the low-middle tropospheric flow is 
southeasterly. During the monsoon period, thunderstorms are an almost daily occurrence in the region. The large 
spatial variability of deep convection over Arizona (Lopez et al. 1997) is directly related to the highly complex 
terrain of the state. Orographic effects tend to produce greater summer precipitation totals in the mountainous areas 
north and east of city (figure 4). In the early afternoon, convective activity peaks on the Mogollon Rim (Watson et 
al. 1994b). What is most interesting about the northeast Phoenix region is the apparent existence of an anomaly 
region (AR) where precipitation totals in the last 50 years appear to be greater than surrounding areas with similar or 
greater topographic relief. 
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ii. Riyadh 


Almost all of Saudi Arabia consists of semi-desert and desert with oases, where half of the total surface is 
uninhabitable desert. Most of the western parts of Saudi Arabia is plateau, the east is lowland, with very hot climate. 
The southwest has mountains as high as 3000 m. Saudi Arabia, one of the driest and hottest countries in the world, 
is roughly located between north latitudes 17 and 31 and east longitudes 37 and 56. Except for the southwestern 
mountains, the average annual precipitation in the Kingdom ranges from 80 mm to 140 mm (Alkolibi, 2002). 
Maximum summer temperatures often exceed 45 degrees Celsius, relative humidity is very low and skies are clear 
most of the time. There is very little precipitation in Saudi's middle (the area known as the Rub Al-Khali desert, or 
the "empty quarter"). But, there is typically more precipitation in the highland and mountainous areas of the 
country. Precipitation increases southwest of the country in the Atlas mountains, as well as to the northeast in the 
Zagros mountains. 

3. Results 

a. Phoenix 

i. Analysis of Pre-Urban and Post-Urban Rainfall in the Phoenix Area 

The pre-urban (1895-1949) and post-urban (1950-2003) precipitation data were grouped such that the mean 
precipitation totals for the monsoon period (July-September) are aggregated. Table 1 summarizes the data for 
selected sites identified in figure 4. Sites were selected to represent the Phoenix region and selected sites to the 
southwest, northwest, northeast, and southeast of the urban area. We analyzed the time series of mean monsoon 
season precipitation for each site to determine if any significant differences between the pre-urban and post-urban 
period were evident. 

Stations located southeast of the Phoenix area like Clifton and Duncan receive a significant amount of the 
precipitation during the monsoon period. Since synoptic-scale monsoon system should cause significantly higher 
precipitation totals in the monsoon zone than in regions around the Phoenix area, comparisons between moderately 
monsoon-impacted stations (i.e. Lower Verde) and strongly monsoon-impacted stations (i.e. “monsoon zone”) 
enable a compelling examination of possible urban-induced increases in summer precipitation totals in the Lower 
Verde or northeast Phoenix region (DB03). The monsoon zone stations are above 1100 m above sea level. The 
Sunflower station, located in the Lower Verde basin northeast of Phoenix is also at significant elevation (> 1 100 m) 
and is the wettest location in our sample. 
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The most interesting result is how the mean precipitation amounts have changed from the pre-urban to post- 
urban periods. Paired t-testing reveal that differences in pre-urban and post-urban precipitation mean totals are 
statistically significant at the 95% level. Regardless of monsoon region or topography, our results identify a 
preferred region of enhanced precipitation in the Lower Verde basin 30-50 km northeast of the Phoenix metropolitan 
area. The only stations in table 1 that have experienced an increase in mean precipitation greater than 10% from the 
pre-urban period to the post-urban period are the Lower Verde Basin stations. The station at Carefree is slightly west 
of Lower Verde but also exhibits a relatively significant increase (9%) over the same period. One sample t-tests 
were conducted for each mean value in the Lower Verde region to see if the differences relative to the other stations 
are statistically significant or occurred by chance. In all tests, the two-tailed P value was less than 0.0001 so by 
conventional hypothesis testing criteria, the difference is considered extremely statistically significant. Though 
DB03 did not investigate the pre-urban period, their results also identified the Lower Verde region as the region of 
elevated median precipitation totals over the period of 1950 to 2000. Stations in other regions have experienced 
essentially no change or slight decreases in mean precipitation totals. Figure 5 is a time series of monsoon period 
(July-September) mean precipitation totals for Sunflower (Lower Verde), Phoenix (Urban), and Clifton (Monsoon). 
The figure suggests that precipitation totals in the Lower Verde station (Sunflower) have trended upward over the 
108-year period as compared to Phoenix and Clifton. 

ii. Analysis of the Relative Magnitude of Precipitation at the Lower Verde Stations 

Another finding of this study is that other stations in the region with comparable topographic relief and similar 
(and in some cases more severe) monsoonal conditions have not experienced significant increases in the post-urban 
period. An analysis similar to DB03 was conducted in which all sites from Table 1 that exhibited a percentage 
increase greater than ten in mean total precipitation were compared with stations in nearby basins and the monsoon 
zone. For consistency with DB03, we analyze the relative differences of the median precipitation totals (Table 2) 
for years 1950-2003 (post urban). 

Results in table 2 are consistent with DB03. It is clear that Sunflower in the Lower Verde is the region’s wettest 
station. In part, this is likely due to Sunflower’s higher elevation although the two monsoon zone stations have 
comparable elevations. Sunflower may also be experiencing a combined urban and topographical effect in the post- 
urban period. Lower Verde’s Horseshoe Dam and Bartlett Dam sites have median precipitation totals greater than 
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sites like Cave Creek and Carefree in the northwest Lower Salt basin even though the elevations are similar (and 
even smaller in the case of the Carefree station). Not surprisingly, the stations in the Monsoon Zone are generally 
wetter than the other sites, except Sunflower, which indicates the influence of synoptic scale monsoon forcing, but it 
is useful to recall from the Table 1 discussion that the Monsoon Zone precipitation totals have remained fairly 
consistent from the pre-urban to post-urban period. 

iii. Hypothesis for the Lower Verde Precipitation Anomaly 

The results of this study were so compelling that we posed the question, “Does the Lower Verde Precipitation 
Anomaly reveal itself under drought conditions? To answer the question, we used the TRMM rainfall product 
(TRP) to construct a mean rainfall rate map for the monsoon period of July-September 2003 for the central Arizona 
region. In figure 6, the TRP estimate of mean rainfall rate reveals an elevated region of rates 30-50 km northeast of 
the Phoenix metropolitan area (approximated by the black rectangle) and the influence of mountainous terrain on 
precipitation as the bluish colors (indicating precipitation) essentially corrolate to the mountains (see figure 4). 
Heavier rainfall rates on the Mexican border likely reflect both topography and stronger monsoonal precipitation. 
These results seem to suggest that the processes leading to the anomaly can exist under the most extreme drought 
conditions. 

Shepherd and Burian (2003) presented evidence that convective anomalies over and downwind of Houston, 
Texas were related to the complex interactions between urban-induced and sea-breeze related dynamics while the 
role of aerosols remain uncertain. Watson et al. (1994) and DB03 have suggested that convective activity in 
Phoenix is enhanced by urban-related processes and the topographic relief. During the late morning to early 
afternoon, thunderstorms are typically confined to the highest terrain (Wallace et al. 1999). Over the central and 
eastern mountains there is a pronounced mid-afternoon maximum of thunderstorm activity directly attributable to 
diurnal heating and mountain-valley circulations. Watson et al. (1994b) have identified an early-evening (i.e. 5 to 7 
pm LST) maximum in deep monsoonal convection in the Phoenix region. 

Phoenix typically exhibits a relatively weak urban heat island (UHI) and sometimes an urban heat sink because 
large amounts of energy are converted into latent heat rather than sensible heat because of the prevalence of irrigated 
lands (DB03). Thielen et al. (2000) showed in numerical simulations that when UHIs are weak, surface sensible 
heat fluxes, convergence, and buoyancy variations that influence rainfall development are most effective at a 
distance from the central heat source. Baik et al. (2001) also showed that when a heat island is present, two regime 


flows are established. One exhibits a stationary gravity wave near the heating region and the other is characterized 
by stationary gravity waves near the heating source and an updraft cell downstream from the heating source. In the 
case of Phoenix, these results suggest that urban-induced convective forcing can be established on the edge of the 
metropolitan area. Figure 7 is a schematic presenting our hypothesized reasons why the Lower Verde anomaly 
exists. First, convection from the mountains is established during the daytime through direct heating and mountain- 
valley thermal circulations. By afternoon, the storms propagate downward and to the west and south releasing 
evaporation-generated density currents or outflow boundaries as hypothesized by Watson et al. (1994b). Figure 8 
uses the TRMM product in a Hovmoller diagram to illustrate the tendency for westward propagation of storms. The 
diagram plots the longitude-time variation of the mean rainfall rate averaged over the 33°-34° latitude box, including 
Phoenix and immediate basin sites. The figure illustrates the tendency for convection to form in the mountains and 
then propagate westward as indicated in this case from August 10-11 2003. Because of the location of the Lower 
Verde region, the outflow boundaries from storms in this region can most effectively interact with the convergence 
and flux forcing on the edge of the city. 

Additionally, figure 3 pointed out that the lower-middle tropospheric flow on active monsoon days is 
southeasterly. The resulting urban-mountain convergent region is continually pumped with monsoon moisture from 
the southeasterly flow. Though speculative, it should be noted that the Phoenix region is likely more moist due to the 
presence of significant and increasing irrigation and agricultural practices. Moisture advection into the preferred 
convergence region by southerly flow regimes could enhance convection. Figure 9 is taken from Gaffen and Ross 
(1999) and illustrates the sign and statistical significance level of station trends in dewpoint temperatures Td for each 
season based on full-day data. The trends are generally positive and statistically significant, particularly at Phoenix 
station. Gaffen and Ross (1999) note that locally, anthropogenic modification of the hydrological cycle is important 
to consider. Within the conterminous United States, the U.S. Geological Survey has estimated that consumptive use 
of water in agricultural irrigation contributes 100 billion gallons of water per day to the atmosphere, compared with 
2800 billion gallons per day from evaporation and transpiration from surface water bodies, land surface, and 
vegetation. In dry regions during the growing season, the ratio of consumptive use to natural evaporative sources 
may be greater, and it is possible that long-term increases in evaporation from irrigated fields may be large enough 
to influence the surface trends at some stations. Numerical simulations using a coupled atmosphere-land surface 
model will be conducted in future studies to test these hypotheses. 
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b. Riyadh 


i. Historical data analysis 

An analysis of GSN mean temperature measurements from 1978 to 1998 for stations at Riyadh (24.4N, 46.4E), 
Madinah (24. 3N, 35.35E), and Al-Jouf (29.79N, 40. IE) indicates that urbanized Riyadh may be experiencing a 
climatological shift towards temperature profile influenced by the well-studied urban heat island. Figure 10 
suggests that Riyadh’s mean annual temperature has trended upward from 1990-1998. It can be argued that the 
temperature at the other two stations have also trended upward over this period but the trend is less clear. These 
results are consistent with Alkolib (2002, hereafter A02) who investigated historical temperature and precipitation 
for the four Saudi weather stations during the period 1961-1997. These four stations were located at Riyadh in the 
central part of the country, Dhahran in the east, Jeddah in the west, and Almadinah in the northwest. A02 employed 
both time series and regression analyses. Five-year centered moving average and regression analyses covering the 
thirty seven-year period 1961-1997 for each of the four stations was applied to the following data: 

A02 found that the Riyadh station does exhibit a noticeable increase in temperature with R2 = 0.39 and F = 
0.001 - a significant value. A20 attributes this increase to the urban heat island associated with rapid growth at 
Riyadh over the 36 year period. To determine whether the noticeable increase in Riyadh’s temperature trend could 
be attributed to a heat island effect, A20 compared temperature records from the old airport to the temperature 
records of the new station at King Khalid International Airport, situated about 20 km northeast of the old airport. 
The stations have similar topographical features. In the analysis, it was shown that while both Almadinah and 
Riyadh show increasing trends, only Riyadh shows a noticeable rise with R2 = 0.35 and F = 0.002 - a significant 
value. This increase is most likely due to a heat island effect. 

Of interest for this study, figure 1 la illustrates the trend in total annual precipitation at Riyadh from the GSN 
data. This result suggests that over the 20-year period, there has been sharp upward trend in rainfall after 1992. 
Total rainfall at Riyadh exhibits an increase from -100 mm in the 1978-1992 period to -300 mm or more after 1992. 
Again, these results are consistent with A20 who found a sharp trend upward in mean annual rainfall at Riyadh after 
1992 (figure 1 lb). In the 1996-1997 period of A20, mean annual precipitation had risen to -180 mm. While other 
stations in the study exhibited slight increases over the same period, there maximum values in the 1996-1997 time 
period are -100 mm (Jeddah), 125 mm, (Dharan), and 95 mm (Almadinah), respectively. This result, though 
inconclusive, suggests that Riyadh has experienced a significant upward trend in precipitation during the same 
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period of time as explosive urban growth such that mean annual (total) rainfall is on the order of 200 mm (450 mm). 
The results are not strongly conclusive in terms of linkages to urbanization, but more research is required to 
understand why there has been a systematic trending upward in precipitation on the Arabian Peninsula. 

ii. Satellite-based analysis 

The dominant precipitation season for Saudi Arabia, relatively speaking, is the period from roughly October to 
April and is primarily convective. During this period, the descending branch of the Hadley cell is not as much of a 
convective suppression mechanism as it is during the summer months. The observational network for precipitation is 
particularly sparse thereby increasing the likelihood that the sparse network will significantly underestimate 
precipitation (Evans et al. 2004). Additionally, Saudi Arabian convective storms appear in a region where relative 
humidity is usually less than 40% and hence a significant proportion of the falling rain evaporates. These factors 
are driving motivations behind the application of satellite-based rainfall estimates in this study. Unfortunately, 
reliable, microwave-based (as opposed to infrared-cloud top techniques) high-spatial resolution (e.g. < 25 km) 
estimates of precipitation from satellites have only become available in the last few years. Therefore, we cannot use 
such data to investigate long-term time series of precipitation. However, we can apply the 3B43 NASA satellite 
product to investigate whether there are any preferred regions for precipitation in and around the central portion of 
Saudi Arabia. 

Figure 12 is the time series of mean accumulated precipitation from the satellite-based product from January 
1998 to May 2003. The mean value is plotted for a 25-km square gridpoint center on Riyadh, Saudi Arabia. The 
figure indicates that the rainy season is during the Northern hemisphere winter. The results also indicate a bimodal 
tendency during this period with corresponding peaks in early and late winter, respectively. It should be emphasized 
that values of rainfall are likely under-estimated because of a “smearing” effect of representing fine convective 
precipitation events on a fairly coarse gridpoint. At a given ground-based rain gauge, values would likely be much 
higher. However, if one can accept the inherent resolution and algorithm flaws with the satellite-based approaches, 
figure 13 indicates the true value of the data. Figure 13a illustrates the ability to resolve a spatial depiction of 
rainfall over a region that is known to be deficient in ground-based rainfall measurements. 

In the figure, it is evident that mountainous terrain west of Riyadh is a primary forcing mechanism for rainfall in 
the region over this period of time. That region has received over 600 mm of rainfall over a roughly 5-year period. 
Of interest for this study, is the relative maximum near 27.75N latitude and 46.75 E longitude. This relative 
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maximum is located roughly 50-75 km north of the urban Riyadh area. Under any type of flow with a southerly 
component, this region would be considered the downwind region of Riyadh. This distance displacement for urban- 
induced rainfall anomalies is very consistent with results from the famous METROMEX study (Changnon et al. 
1981) and more recent studies by Shepherd et al. (2002). it is tempting to attribute the anomaly to topography but 
figure 13b suggests that this particular regions is does not have significant relief. Recalling a previous discussion, 
Thielen et al. (2000) and Baik et al. (2001) have shown that weak or small heat islands result in surface sensible heat 
fluxes, convergence, and buoyancy variations that influence rainfall development being most effective at a distance 
from the central heat source. According to Sahin and Aksakal (1999) prevailing low-level wind flow during the 
rainy season is North-Northwesterly. Under this scenario, it is possible that the low-level prevailing flow interacts 
with the weak urban-induced circulations on the northern fringe of Riyadh to establish a preferred convergence or 
rising branch region. Additionally, irrigation has likely altered the moisture availability in ways similar Phoenix 
(e.g., figure 9). However, at this point, this hypothesis has not been confirmed with numerical simulations. 

4. Discussion and Conclusions 

This study has employed a unique 108-year precipitation data record, global observing network and satellite- 
based to identify statistically significant anomalies in rainfall downwind of the Phoenix, Arizona urban region and a 
possible urban-induced anomaly just north of Riyadh, Saudi Arabia.. Our analysis showed that during the monsoon 
rainfall season (July, August, September) locations in the northeastern suburbs and exurbs of the Phoenix 
metropolitan urban area (the so-called “anomaly region”), primarily the Lower Verde Basin, have experienced 
statistically significant increases in mean precipitation of 12 to 14 percent from a pre-urban (1895-1949) to post- 
urban (1950-2003) period. Additionally, mean and median post-urban precipitation totals in the anomaly region 
were significantly greater, in the statistical sense, than regions west of the city and in nearby mountainous regions of 
similar or greater topography. Further analysis of satellite-based rainfall totals for the summer of 2003 also revealed 
the existence of the anomaly region even though the arid southwest experienced severe drought during the period. 
More significant orographic features do not correspond to the location of the anomaly region, and statistically 
significant precipitation changes are found in the anomaly region when pre- and post- urban data are analyzed. For 
the Saudi Arabian case, results from our analysis of gauge data did not inclusive link any precipitation trends with 
urbanization. However, satellite-based analyses revealed a relative maxima in accumulated rainfall 50-100 km north 
of Riyadh. 
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Future work will utilize coupled atmosphere and land surface models to test the hypothesis that urban- 
topographic dynamics and increased moisture from irrigation practices alter precipitation in arid regimes. The role 
of aerosols must also be investigated. The implications of this research are wide. The results continue to suggest 
that weather and climate models must adequately resolve the urban environment (e.g. land use, aerosols). The 
results also indicate that even a fairly significant large synoptic-scale process (e.g. monsoon) might possibly interact 
with local and regional mesoscale circulations. Finally, our results highlight how sensitive the water cycle can be to 
anthropogenic forcings, even under extreme drought or arid conditions. 
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Table 1-Comparison of pre-urban (1895-1949) and post-urban (1950-2003) mean precipitation 
totals in the Phoenix and surrounding basin areas. 


Location 

Latitude 

Longitude 

Region and 
Elevation (m) 

Pre-Urban ( 1 895- 
1949) Mean 

Precipitation (mm) 

Post-Urban 

(1950-2003) 

Mean 

Precipitation 

(mm) 

Percent 

Change 

Phoenix 

(downtown) 

33.4 

111.9 

Urban 
(610 m) 

23.66 

23.39 

-1.1% 

Phoenix 

(SW) 

33.0 

112.5 

Urban 
(304 m) 

22.74 

22.38 

-1.5% 

El Mirage 

33.6 

112.3 

Urban 
(—610 m) 

21.00 

21.63 

3.0% 

Buckeye 

33.3 

112.5 

Urban 
(304 m) 

23.26 

22.28 

-4.2% 

Duncan 

32.7 

109.1 

Monsoon 
(1160 m) 

41.41 

43.04 

3.9% 

Clifton 

33.0 

109.2 

Monsoon 
(1150 m) 

46.04 

46.20 

0.3% 

Kelvin 

33.1 

110.9 

Middle Gila 
(550 m) 

37.73 

37.15 

-1.5% 

Florence 

33.0 

111.3 

Middle Gila 
(460 m) 

26.70 

28.05 

3.2% 

Stewart 

Mtn 

33.5 

111.5 

East Lower 
Salt (430 m) 

33.68 

34.78 

3.2% 

Roosevelt 

Dam 

33.6 

111.1 

East Lower 
Salt (650 m) 

38.51 

40.93 

6.2% 

Sunflower 

33.9 

111.4 

Lower Verde 
(1160 m) 

58.95 

65.91 

11.8% 

Horseshoe 

Dam 

33.9 

111.7 

Lower Verde 
(615 m) 

34.21 

38.19 

11.6% 

Bartlett 

Dam 

33.8 

111.6 

Lower Verde 
(500 m) 

29.48 

33.59 

13.9% 

Cave Creek 

33.7 

112.0 

NW Lower 
Salt (490 m) 

21.64 

22.91 

5.8% 

Carefree 

33.8 

111.9 

N W Lower 
Salt (760 m) 

27.55 

30.05 

9.0% 
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Table 2-Relative Differences (mm) between Median Summer Precipitation Totals at Lower Verde 
Stations and other study station, and (Years 1950-2003). 


Location 

Bartlett Dam 
(500 m) 

Horseshoe Dam 
(615 m) 

Sunflower 
(1,160 m) 

Region 

Phoenix (downtown) 
(610m) 

+7.05 

+12.97 

+40.03 

urban 

Phoenix (SW) (304 m) 

+8.41 

+14.33 

+41.39 

urban 

El Mirage (~610 m) 

+11.64 

+17.56 

+44.62 

urban 

Buckeye (304 m) 

+10.68 

+16.60 

+43.66 

urban 

Duncan ( 1 , 1 60 m) 

-11.43 

-5.51 

+21.55 

monsoon 

zone 

Clifton (1,150 m) 

-19.91 

-12.12 

+14.94 

monsoon 

zone 

Kelvin (550 m) 

-7.06 

-1.14 

+25.92 

middle gila 

Florence (460 m) 

+3.45 

+9.37 

+36.43 

middle gila 

Stewart Mtn (430 m) 

-2.73 

+3.19 

+30.25 

lower salt 

Roosevelt Dam (650 
m) 

-9.9 

-3.98 

+23.08 

lower salt 

Cave Creek (490 m) 

+7.83 

+14.25 

+41.31 

nw lower 
salt 

Carefree (760 m) 

+1.8 

+7.79 

+34.85 

nw lower 
salt 
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Figure 1-LANDSAT images illustrating the growth of the Phoenix metropolitan area from 1973 to 1989. 
The yellow outline encompasses the urban region. Red colors are crops or vegetation. 



Figure 2-LANDSAT and ASTER satellite images of the growth of Riyadh. Saudi Arabia from 1972- 
2000. Grey areas represent urban surfaces, and red areas represent irrigated land. 
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Figure 5-Location of stations examined relative to central Arizona basins. Urban and irrigated cropland 
were identified using a 1999 Landsat Enhanced Thematic Mapper scene (figure taken from Diem and 
Brown 2003) 
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Figure 6-Mean rainfall rates for the 2003 monsoon period (July-September) derived from the real-time 
TRMM 3B42 rainfall product. Black rectangle is the approximate location of the Phoenix urban area. 
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Figure 7-Schematic of the hypothesized urban-topographic interactions in the Phoenix region during a 
typical monsoon season. 
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Figure 8-TRMM 3B42RT Hovmoller diagram illustrating the tendency for Phoenix area monsoon 
convection to develop in the mountains east of Phoenix then propagate down in the Phoenix-urban valley 
area with a westward component. The plot represents mean rainfall rate for a convective period of 09 
UTC 10 August to 09 UTC 1 1 August. 
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Figure 9- Sign and statistical significance of trends, b. in dewpoint (K decade- 1 ) for (a) winter, (b) 
spring, (c) summer, and (d) fall based on full-day data. Statistical significance, indicated by the two-sided 
probability that the rank-order correlation between dewpoint and time is not different from zero, is shown 
for the 0.05 and 0.01 confidence levels (from Gaffen and Ross 1999). 
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Figure 10-GCOS Surface Network temperature records (1978 to 1998) for Madioah. Al-jouf, and Riyadh. 
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Figure 11 a. (top) GCOS Surface Network record of total annual precipitation at Riyadh (1978 to 1998 ) 
b. (bottom) Trend of mean annual precipitation at Riyadh from 1961 to 1997 (following Alkolib (2002)). 
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Figure 13 a. (top) 3B43 Satellite-based estimate of mean precipitation (mm) from January 1998 to May 
2003. The black circle is the approximate location of the Riyadh urban area. b. (bottom) Relief map 
centered on Riyadh with significant topographic relief to the west of Riyadh (map courtesy of ESRI, 
EarthSat, and National Geographic). 
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Popular Summary 


Evidence of Urban-Induced Precipitation Variability in 
Arid Climate Regimes 

Submitted to Journal of Arid Environments 

J. Marshall Shepherd 
NASA/GSFC 

Water is essential to life in the Earth system. The water cycle components that sustain life are 
becoming more scarce and polluted. The most recent (1999-2004) drought experienced in the 
southwestern United States is the seventh worst in the approximately 500-year proxy tree-ring 
record. As a result, many regions contemplated “drought emergencies” in which severe water 
restrictions are implemented. Though larger weather and climate processes likely control drought 
processes, there is increasing evidence that anthropogenic or “human-related” activities can 
significantly alter precipitation processes. Urbanization is an example of anthropogenic forcing. 
Recent studies continue to provide evidence that urban environments can modify or induce 
precipitation under a specific set of conditions. 

Arid and semi-arid regions of the southwestern United States and other parts of the world are 
rapidly developing and placing greater demands on the environmental system. In the past fifty 
years, Phoenix has expanded from a predominantly agricultural center to an urbanized region 
with extent 700 percent larger than its size in the middle of the twentieth century. Riyadh’s 
population grew from about a half million people in 1972 to almost two million by 2000. Saudi 
Arabia experienced urbanization later than many other countries; in the early 1970s its urban- 
rural ratio was still about 1:3. By 1990 the ratio had reversed to about 3:1. In the mid-1970s 
Riyadh’s population was increasing by about 10 percent a year. Irrigation also significantly 
increased between 1972 and 1990 southeast of Riyadh. 

The study employs a 108-year precipitation historical data record, global climate observing 
network observations and satellite data to identify possible anomalies in rainfall in and around 
two major arid urban areas, Phoenix, Arizona and Riyadh, Saudi Arabia. It provides statistically 
sound evidence that rainfall distribution and magnitude is statistically different in post-urban 
Phoenix than in pre-urban (1900-1950) Phoenix. The study hypothesis that a complex 
interaction between the city landscape, irrigated lands, and nearby mountains have created 
preferred regions for rainfall development. The study also provides early evidence that rapidly 
urbanizing parts of the arid Middle East may also be experiencing different precipitation regimes 
in response to urbanization and irrigation. 


